Bogoliubov angle and visualization of particle-hole mixture in superconductors. 
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Superconducting excitations — Bogoliubov quasiparticles — are the quantum mechanical mixture 
of negatively charged electron (-e) and positively charged hole (+e). Depending on the applied volt- 
age bias in STM one can sample the particle and hole content of such a superconducting excitation. 
Recent Scanning Tunneling Microscope (STM) experiments offer a unique insight into the inner 
workings of the superconducting state of superconductors. We propose a new observable quantity 
for STM studies that is the manifestation of the particle-hole dualism of the quasiparticles. We call 
it a Bogoliubov angle. This angle measures the relative weight of particle and hole amplitude in the 
superconducting (Bogoliubov) quasiparticle. We argue that this quantity can be measured locally 
by comparing the ratio of tunneling currents at positive and negative biases. This Bogoliubov an- 
gle allows one to measure directly the energy and position dependent particle-hole admixture and 
therefore visualize robustness of superconducting state locally. It may also allow one to measure the 
particle-hole admixture of excitations in normal state above critical temperature and thus may be 
used to measure superconducting correlations in pseudogap state. 

PACS numbers: Pacs Numbers: 



I. INTRODUCTION 

The dual particle-wave character of microscopic ob- 
jects is one of the most striking phenomena in nature. 
This dualism is ubiquitous in the microworld. Most no- 
tably, the two-slit experiments of Stern and Gerlach re- 
vealed the interference and, hence, the wave nature of 
electrons. In the condensed matter systems, such explicit 
visualization of the wave nature of the constituent elec- 
trons was missing until just recently. The breakthrough 
came when the researchers from the IBM labs realized 
that the best way to elucidate the electrons inside a ma- 
terial is to place an impurity in an otherwise perfect crys- 
tal structure. By building corrals of the impurities on 
the clean surface, and observing the generated patters 
through the scanning tunneling microscope (STM), the 
experimenters were able to demonstrate the laws of the 
wave optics using the conduction electron wavesji^^i^ 

The analog of the conduction electrons in the supercon- 
ductors are the quasiparticles. Unlike electrons, the su- 
perconducting quasiparticles do not carry definite charge. 
The same quantum mechanical dualism is at play when 
one considers the Bogoliubov quasiparticles in supercon- 
ducting state: the quasiparticle is a coherent combina- 
tion of an electron and its absence ("hole"). Particle- 
hole dualism of quasiparticles is responsible for a variety 
of profound phenomena in superconducting state such as 
Andreev reflection, the particle-hole conversion process 



that is only possible in superconductor. 

In this paper we propose a technique to reveal this co- 
herent particle-hole mixture locally. In order to discuss 
the particle-hole mixture we introduce a quantity that 
parametrizes the mixture in terms of an angle, we call 
this angle a Bogoliubov angle (BA), see Fig. [TJ We ar- 
gue that STM measurements allow one to visualize the 
Bogoliubov angle maps and thus to reveal particle hole 
dualism. Bogoliubov angle maps as a function of po- 
sition and energy offer a tool to investigate strength of 
superconducting state locally. 

Bogoliubov showed that in order to obtain natural ex- 
citations in the superconducting state one needs to use a 
linear combination of particle in hole excitations with the 
coherence factors Mn(ri) and i'„(rj). They describe the 
unitary transformation from particle and hole operators 
to quasiparticles that are: 



(1) 



with the constraints that /(ir(|u„(ri)p - 
for any n (normalization) and X^n (I ) P + 1 ) P ) — 
1 for any i (orthonormality), i being the site index on our 
lattice. 

In the normal state, either u„(r.i) or Vni^i) are iden- 
tically zero and there is no mixing between the particle- 
and hole-component of Bogoliubov quasiparticle. Once 
superconductivity sets in, the mixing between these com- 
ponents develops. This mixing strength can be repre- 
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sented by 

e„(rO=arctan((^^4^)'^') (2) 

which is a central quantity we are interested in. We define 
this quantity as a Bogoliubov angle. The high resolution 
STM allows us to study the spatial dependence of the 
BA for the states whose energy can be selected by tuning 
STM bias. 

hole 




FIG. 1: Circle parametrizing Bogoliubov admixture angle is 
shown. For O = 0, 7r/2 the mixture reduces to purely hole-like 
and particle-like state. At arbitrary angle one deals with true 
Bogoliubov quasiparticles. 

Note, we intentionally do not simplify the expression 
in Eq.© for the reasons that will be clear in the next 
section. It represents a local mixture between particle 
and hole excitations for an eigenstate n at a given site 
i. For example, for 8„(ri) = the Bogoliubov excita- 
tion will be a hole. In the opposite case of 8„(ri) = 7r/2 
quasiparticle is essentially an electron. The angle that 
corresponds to the strongest admixture between particle 
and hole is 8„(ri) = 7r/4 = 45°. Obviously, in case of 
inhomogeneous state the BA is a function of a position 
where it is measured and also is a function of energy E. 
We suggest a way to visualize the BA maps that allow us 
to develop a more detailed understanding of the super- 
conducting state. Previoisly the alternation of coherence 
factors u, v as functions of position near impurities has 
been discusses Here we expand this discussion 

by introducing BA. We also will focus more on the spon- 
taneous inhomogeneity and not the impurity states that 
were the focus of previous studies. 

The ideas presented here are quite general and are ap- 
plicable to a variety of superconductors, including con- 
ventional superconductors. Imaging of BA can be per- 
formed in any inhomogeneous state. One can investi- 
gate BA in a variety of states, including vortex state and 
normal state with superconducting correlations, e.g. so 
called pseudogap (PG) state^i^. To illustrate this ap- 
proach we will use the local STM data obtained on high- 
Tc superconductor, namely on Bi2Sr2CaCu20a+5 mate- 
rial. 



The plan of the paper is as follows. We first present a 
general theoretical background and define BA from the 
local tunnehng conductance measurements dI/dV{r,V) 
at different bias values V. Then we describe the numer- 
ical results for the calculation of BA. 



II. THEORETICAL DISCUSSION 

To illustrate the point about BA, we can look at the 
uniform BCS case first. Using Bogoliubov quasiparticles 
one can introduce BA as: 

0^ = -ctan[(M^)V^], (3) 

with the conventional coherence factors, see Fig. [21 

We now turn to inhomogeneous problem where we can 
no longer use translational invariance and plane waves as 
a basis. It is well known that the wave function for a BCS 
superconductor is well captured by a following mean field 
wave function: 

n 

where n is the eingenvalue index with respective single 
electron creation operators cj^ ^ in the state n with spin 
s. Real space description used here is necessary in case 
we consider the effects of disorder in single particle po- 
tential. Either kind of disorder breaks translational sym- 
metry and the real space representation is more natural 
in this case. The pair wave function is captured in the 
relative position dependence of u„(ri — ri'),i;„(ri — r^/) 
components of Bogoliubov spinor. The self-consistently 
defined gap amplitude 

A(ri,rj) = Vint^Un{ri)vnirj){l - f{En)), (5) 

n 

with Vint being interaction, and f{En) being the Fermi 
distribution function for a given quasiparticle excitation 
spectrum En- Here En are defined to be positive only. 
We assume that pairing interaction couples only nearest 
neighbors on the lattice Yj , . 

For a next step in our discussion it is necessary to 
introduce a tunneling conductance as measured by local 
STM tunneling. Introduce the tunneling conductance on 
positive and negative bias E — ±e|y| as: 

dI/dV+in,E) = Fiz,\eV\)Y,K,s\Hn)i-l)f'iE - En), 

n,s 

dI/dV-{v,,E) - F{z,-\eV\)Y,\vn,s\'{y^){-l)f'{E + En), 

n,s 



where i^(z, ±|ey|) is a function that measures the ma- 
trix elements for tunneling as a function of voltage bias 
and tip distance z. Hereafter we assume that it is a 
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smooth function of energy and at small energies E ~ 
10 — lOQmeV it is a constant. f{E) is Fermi distribu- 
tion function. At very low temperatures f'{E) becomes 
a nearly 5{E) function, a fact that we will use often. We 
can simplify the formulas if we introduce the density of 
states (DOS) for quasiparticles: p{E) = X)„ s 5{E~ En). 
For simplicity we will assume particle-hole symmetry in 
the normal state. 

Then, for a given eigenspectrum and eigenfunctions 
Uni^i — ri/),w„(ri — r^') we can rewrite Eq.® as: 



dI/dV+{r,,\eV\) = - J dEp{E)\uE^. 
E>0, 

dI/dV^ir,,\eV\) = - [ dEpiE)\vE,s\\r,)f{\eV\+E), 



\r,)f'i\eV\~E), 



E <0. 

(7) 

Hence the ratio of dl/dV, taken at the same \E\, that we 
label as Z{ri, \eV\), will be 



Z{r„\eV\^E) 



dI/dV+{r„\eV\) 
dI/dV-{n,\eV\) 



\uE,s\'^iri) 

\vE..s\Hr^) 



= tan2e(r„S), (8) 



where the last step is taken assuming that there are few, 
often one state, that contributes to the summation in 
Eqs.® a energy E — \eV\. Then Eq.® can be inverted 



Q{ri,E) = arctan[( 



dI/dV+{r,,\eV\) 
dI/dV-{r„\eV\)' ^ 



(9) 



this result along with Eq.([2|) is the main result of this 
section. It allows a direct determination of Bogoliubov 
angle 6(ri, E) from the experimentally measured tunnel- 
ing conductances at positive and negative bias. 

BA as a measure of particle-hole admixture appears 
naturally in the Anderson mapping^'' of BCS model on 
the effective spin model. We briefly recall the mapping 
in Appendix A. 

To visualize the local quasiparticle states we employ 
the Scanning Tunneling Miscroscopy (STM) technique. 
Crucial aspect of the electron tunneling into the su- 
perconducting state that makes it qualitatively differ- 
ent from the tunneling in conventional metals is that 
the STM tip contains only the regular electrons which 
carry a unit of charge (-e). We can inject either elec- 
trons or holes in superconductor. On the other hand 
as was pointed out early on starting with Bogoliubov, 
quasiparticles that live inside the superconductor do not 
possess a well-defined charge. Upon entering the super- 
conductor, an electron/hole that arrived from the normal 
STM tip must undergo a transformation into the Bo- 
goliubov quasiparticles native to the superconductor—. 
Hence electrons that are injected or extracted form su- 
perconductor would need to be "assembled" from Bo- 
goliubov excitations. At any site and at specific bias this 



conversion into particles and holes will depend on relative 
weights Un(ji), Vn{vi). Hence the intensity of a tunneling 
signal will depend on these coherence factors. 

Qualitatively, the spatial distribution of tunneling in- 
tensity can be understood as follows. Respective ampli- 
tudes of particle and hole parts of the Bogoliubov quasi- 
particle, are Uniji) and Vn(vi) for site i and for particular 
eigenstate n. Consider now a site where, say, M„(ri) is 
large and close to 1. It follows therefore that for the same 
site the Vni^i) would have to be small, since the nor- 
malization condition is almost fulfilled by |un(ri)P term 
alone. Similarly, for the sites where w„(ri) has large mag- 
nitude, u„(rj) would have to be small. Recall now that 
large Un(vi) component would mean that quasiparticle 
has a large electron component on this site. Hence the 
electron will have large probability to tunnel into super- 
conductor on this site and the tunneling intensity for elec- 
trons (positive bias) will be large. Conversely, for those 
sites the hole amplitude is small |f (ri)| <C |M(ri)l ^^'^ the 
hole intensity (negative bias) will be small. Similarly, 
for sites with large hole amplitudes |w(ri)| » |u(ri)| the 
electron amplitude will be suppressed and this site will 
be bright on the hole bias. We observe alternation of the 
form: 



|x;„(r,)|2 ~ 1, |M„(ri)|2 < 1, 
|w„(r,)p < 1, |M„(ri)p ~ 1. 



(10) 



Therefore if there is a particular pattern for the large par- 
ticle amplitude (sampled on positive bias) on certain sites 
i, the complimentary pattern of bright sites for hole tun- 
neling (on negative bias) will develop as a consequence 
of the inherent particle-hole mixture in superconductor. 
This antiphase behavior is a clear indication of the "nat- 
ural quasiparticles" having both particle and hole char- 
acter. It is the main effect that can be visualized by 
considering Q(r i,E) maps. Antiphase shift in positive 
and negative bias intensity is ubiquitously seen in tun- 
neling spectra. The "antiphase" behavior of the compo- 
nents |it£;(ri)P' is explained here as a case of 
BA changing from particle to hole-like configuration on 
alternating sites. We see that this is the case in our nu- 
merical simulations, (see Numerical Simulations below), 
without any need to assume that only one state domi- 
nates the sum over states in Eq.(l7]). So the phenomenon 
is more general. We find it easiest to explain assuming 
only one term dominating. But given numerical results 
it holds for broader cases. 

We discuss it in more details below when we turn to 
Q{r i,E) maps. 



A. Particle-Hole asymmetry of normal state 

The question of the underlying band particle-hole 
asymmetry often comes up in these materials at low dop- 
ing. One way to " factor out" this asymmetry that is ex- 
trinsic to the particle-hole mixture measure, is to factor 
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out the normal state conductances; namely one can take 
a ratio of dI/dV±{ri,V,T) to their proper normal state 
values at high temperatures T > Tc/. 



dI/dV±{r,,E, T) 



dI/dV±{r,,E, T) 
dI/dV±{ri,E,T>T,) 



(11) 



This procedure will factor out the particle hole asymme- 
try for the underlying band and will allow more direct 
measure of particle-hole asymmetry. 



III. IMAGING BOGOLIUBOV ANGLE IN 
NORMAL AND PSEUDOGAP STATE 

The BA, as defined is not sensitive to the SC quan- 
tum phase fluctuations. Indeed BA is defined as a 
function of ratio of the |u_E(i"i)P/|i'_E(ri)P- Therefore 
Q{ri,E) can be defined even in the presence of such 
phase fluctuations^^. Thus we propose that the discus- 
sion about BA may be extended to the normal state. 

Imagine we are approaching a normal state of super- 
conductor by warming it up. We can see that there will 
be temperature dependence of the BA. There is no rea- 
son to expect an abrupt termination of SC correlations 
as one crosses Tc. Remnant superconducting correlations 
are present above Tj^^ and hence one can still have ex- 
citations that will have a particle-hole admixture. The 
difference will be that we are no longer in the state with 
well defined superconducting Josephson phase. 

To illustrate this point consider Bogoliubov-Valatin 
transformation in the presence of phase fluctuations: 

7n,T(i"i) = "«(ri)c„^| -I- exp(i0(ri))w„(ri)cjj^^. (12) 

We then can use the same definition for BA, Eq.([21) in 
this case even in the presence of random Josephson phase 
(p. Since one uses amplitudes of u, v, spatial phase disor- 
der does not enter into Q{ri,E). So for the frozen and 
presumably for the slowly varying in time phase fluctua- 
tions once can use the BA as defined and image the local 
particle-hole admixture in the normal state. 

One would need to take care of thermal broadening 
of the tunneling characteristics at higher temperatures. 
Namely one could divide the tunneling characteristics by 
derivatives of the Fermi thermal distribution function : 



&PGiri,V) = arctan[[ 



dI/dV+{r,,V)f{E + \eV\) 
dI/dV-ir,,V)f'iE-\eV\y ^ 



(13) 



The problem of dynamic phase fluctuations in the state 
with superconducting fluctuations is complicated. More 
detailed analysis would require a specific model for the 
dynamics of the superconducting phase. An approach 
to phase fluctuations in PG state using localized Cooper 
pairs state with no long range phase coherence was ad- 
vocated in-^, using STM dataS. 



One can also study the behavior BA for other states, 
such as flux phaseiS. state and density wave statesii. 
Consider density wave states, e.g. d- density wave state 
(DDW). DDW is often mentioned as a possible state that 
can explain PG^^. In any density wave state, including 
DDW, particle-hole symmetry is violated and the poles 
of single particle excitations are not appearing in pairs 
symmetrically around chemical potential. Therefore sin- 
gle electron tunneling DOS does not have the components 
that appear symmetrically at positive and negative bias. 
If there is a particle hole symmetric spectrum for DDW 
state it can occur only as a special case at one doping 
level. 

Absence of particle-hole symmetry will be easily de- 
tected by BA as it will tend to pure hole or particle 
angle, 8 — > 0,7r/2. Thus we think BA can be used as 
a spectroscopy tool to detect presence/absence of super- 
conducting correlations in normal state. Another inter- 
esting question to address is how BA behaves upon rising 
temperature. At low energy it will be close to 7r/2 but 
then it can quickly move away to indicate purely particle 
or hole states at T > Tc for non-pairing PG state. 

These questions go beyond the scope of this paper and 
will be addressed in separate publication. 



IV. EXPERIMENT 

In order to visualize the BA, we have per- 
formed an experimental investigation of the Spectro- 
scopic Imaging Scanning Tunneling Microscopy (SI- 
STM) measurement on high temperature superconduc- 
tor Bi2Sr2CaCu208+*— . A single crystal of Bi-2212 
grown by floating zone method, is hole-doped by intro- 
ducing non-stoichiometric oxygen atoms per unit cell, 
and its hole concentration is adjusted for slightly over- 
doping {Tc =89K). The crystal is cleaved in the ultra- 
high vacuum and immediately inserted into the STM 
head at T —A.2K. To show the BA, we acquired local 
density of states (LDOS) images by measuring the STM 
tip-sample differential tunneling conductance g{f, V) = 
dl/dV at each location T^and bias voltage V. Since 
LDOS(r, E — eV) (x g{f, V), energy and position depen- 
dence of the LDOS is obtained. 

In Fig. [3^, we show 54nm g{f, V) map at y=-16mV 
measured on the Bi-2212 surface, showing the spatial 
modulations which is interpreted as an interference of 
the Bogoliubov quasiparticlesi^. A Fourier transform 
of g{f,—18mV) in the inset of Fig. [3] exhibits several 
fourier spots corresponding to the period of modulation 
in real space. These observations are consistent with pre- 
vious report a^°i^^ . Although similar modulations are vis- 
ible in g(r, -|-18mV) shown in Fig. [3] which is the same 
FOV (field-of-view) as Fig. [3^, one can notice that the 
spatial phase of these modulations is different. dl/dV- 
spectra which are averaged over the regions with the 
same gap size, where g{f, — 18mV) in Fig. [3^ is slightly 
higher/lower then the average value, see black/red curves 
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in Fig. [S]:. Overall feature of the spectra taken with dif- 
ferent intensity of (7(r, — 18mV) with the same gap are 
almost identical, however, the significant differences at 
low energies in the spectra are seen in the Fig. It is 
obvious that the spectrum with relatively higher ampli- 
tude at the negative sample bias has the relatively smaller 
amplitude at the positive sample bias (and vice versa). 
This implies that the particle and the hole in the super- 
conducting state are entangled each other. 

In Fig. 2^, we calculate a local BA 9(ri, V) by taking 
ratio of the positive and the negative sample bias g{f, V), 
using the following simple formulas, 

e(r„y) = arctan(\/Z). (15) 

Taking ratio has an advantage to cancel out the un- 
known matrix element involved in dl /dV^^ . Fig. 3^ is 
the BA map at 1^ = 18mV with its fourier transform, 
and we found that the 9(ri,18mV) shows spatial mod- 
ulations as well as dl/dV map in the Fig. [31 but with 
more stronger contrast. As seen in the Fig. [3)1, ampli- 
tude of dl/dV between positive and negative bias are 
anti-correlated, so that the taking ratio enhances such 
structure, namely, spatial modulations. BA map is es- 
sentially different from the dl/dV map, since BA map 
exhibits the degree of spatial particle-hole mixture of the 
Bogoliubov quasiparticles. However, as evidenced by the 
fourier transform of G(ri, ISmV) ( Fig. [Dd), fourier pat- 
tern is qualitatively the same as those of Fig. [5^ and b, 
indicating that the period of the existing modulation in 
the BA map is similar to dl/dV modulations. This simi- 
larity supports the claim that the local particle and hole 
amplitude are modulated by scattering. Taking the ratio 
of dl/dV in Eq.(|l4l) and taking the BA map in Eq. fTS)) 
can therefore be important new tools to search for the 
true spatial modulations and individual fourier spots in 
the electron density of states. 

In Fig. [It, we show the distribution of BA at energy, 

V = 18mV which is peaked at = 43°, not exactly at 
45°. One possibility is that the apparent shift of the dis- 
tribution is caused by asymmetric background in the tun- 
neling spectrum^-' that is sampled more at higher voltage 
as we shall see in the Fig. [5]f. To visualize the particle- 
and hole-like regions more clearly, line cuts of BA at 

V =18mV as well as dl/dV atV ^ +18mV and -18mV, 
along the trajectry shown in Fig. [3^, are exhibited in 
Fig. [4jl. For simplicity, wc only fucus on the specific q- 
vector in Fig. [JJa, so that the line profiles are taken from 
the fourier-filtered dl/dV and 8(ri, V = 18mV) with 177- 
vector highlighted by red circle in Fig. [4j3. Particle- and 
hole-like regions are spatially modulating along the line 
and clearly show the anti-phase behavior in modulation 
between dl/dV at -|-18mV and -18mV. 

Figures [HI a-e) show the Q{ri,V) maps for various 
bias voltages and their fourier transforms. With in- 
creasing energy, the periods of modulation in real space 



change, and corresponding fourier spots in the inset of the 
Fig. [5Ka-e) move, following the octet modeliS,, and these 
observations are consistent with previous reports22i2i. In 
addition to the period of modulation, one can immedi- 
ately notice that the pattern of the spatial modulation 
changes. At low energies (Fig. [S^ and b), spatial modu- 
lations are visible all over the field of view. On the other 
hand, at V =34mV (or V >34mV), such modulations 
tend to be visible in the restricted area. This difference 
implies that the different type of scattering might kick in 
at V =34mV (or V >34mV). 

In Fig. [SF, we show the 2D distribution of the BA 
in which distributions are normalized at each energies. 
The spatial change in the BA map seems to occur as a 
crossover, and it can be realized by deviation of the BA 
from 8 = 45° in the Fig. [Sf. The energy which differ- 
entiates the spatially coherent excitations and the local- 
ized excitations is estimated as ^26mV (less than mean 
A ^40mV) where BA starts to monotonically decrease. 

The visualization of the BA shed light to under- 
stand the quasiparticle excitations in the superconduct- 
ing state. The interferences of the Bogoliubov quasipar- 
ticles can be understood as a spatial variation of relative 
weight of the particle and the hole amplitude which is 
represented by the BA. The BA can be a measure of 
the energy scale of the coherent excitations which split 
the type of the modulation structure in real space. And, 
since the spatial modulations in the electronic structure 
are revealed much more clearly in the BA map, this pro- 
vides an excellent new technique to determine the mo- 
mentum space (g-space) electronic structure using SI- 
STM. Hanaguri et al. have recently demonstrated the 
power of this technique with the discovery of the inter- 
ference of the Bogoliubov quasiparticles in Na-CCOCSi. 

V. NUMERICAL SIMULATIONS 

We implement simple but realistic model of the opti- 
mally doped cuprate superconductor with disorder. We 
use a simple BCS solution to illustrate the approach on 
how one can visualize the supercondcuting admixture of 
particles and holes in the natural Bogoliubov excitations 
in superconducting state. Even though the model is sim- 
plistic the approach itself is quite general. 

To model the high-temperature superconductors we 
utilize the highly-anisotropic structure of the cuprates 
and focus on a single layer of the material. In the simpli- 
fied model, the conduction electrons live on the copper 
sites, i, and can hop to the neighboring sites, j, with a 
certain probability measured by the quantity t. In addi- 
tion to that, the electrons that occupy the neighboring 
sites feel mutual attraction of a strength Vint ■ Formally, 
this model is represented by the Hamiltonian, 

Hq = -t ^ cl^Cja-Vint ^ niUj, (16) 
<i,j>,a <i;j> 

were a quantum-mechanical operator creates an elec- 
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tron on site i, the operator Cja- ehminates an electron 
from the site j, and rii = c]^|Ci| + cl^Cii represents the 
electron density on site i. The electron spin, a, can point 
up or down. This model, referred to as the t — V model, 
is known to produce the d-wave pairing for the electron 
densities close to one electron per lattice site, and has 
been successfully used to describe strong impurities in 
a d-wave superconductor^^. The local impurity is intro- 
duced by modifying the electron energy on a particular 
site. The corresponding correction to the Hamiltonian is 

H,^p^V"^P{n,T+n^^). (17) 

This term is the potential part of the impurity energy 
that couples to the total electronic density on site i. We 
solve the impurity problem in the Hartree-Fock approx- 
imation, which replaces the two-body interaction in Hq 
with an effective singe-electron potential. Our goal is to 
use Vimp to investigate the spatial distribution of BA as 
a function of position and energy. 

Using the Bogoliubov-Valatin transformation to the 
quasi-particles operators 7„o-, 

n 
n 

and the mean-field approximation, one can diagonalize 
the Hamiltonian Eq. (jl6p . The quasiparticle amplitudes 
on lattice sites (uniri) , Vn{ri)) have to satisfy inhomoge- 
neous Bogoliubov-de Gennes equationsS^: 




where the kinetic operator ^ and superconducting order 
parameter A can be represented as: 

s 

Avnir,) = J2^s{n)vn{r, + d), (20) 

5 

where S = ±x, ±y are nearest neighbor vectors for a 
square lattice. 

We solve Eq. fTO)) together with the self-consistency 
condition: 

Aa(rO = ^EK(^'' + ^)<('-^) + 

n 

u„(r,)<(r, + S)) tanh(^„/2fcBT), (21) 

where the summation is over the positive eigenvalues En 
only. 

For a square lattice system with Lx L lattice sites, the 
solution of the Bogoliubov-de Genes equations Eq. ([20]) is 
equivalent to the eigenproblem for a 2L^ x 2L^ matrix. 



In order to minimize the boundary effects we assume pe- 
riodic boundary conditions in both x and y directions. 

We have performed numerical simulations on a square 
lattice 32 X 32 at T = 0. We assume 40 impurities ran- 
domly placed on the lattice, each impurity has strength 
yimp _ corresponds to approximately 3% doping. 

We set Vint = and a half-filled band /i = 0. 

The results of our numerical simulations are summa- 
rized in Figs. I6|7|81 where we show three panels of plots 
for calculated local tunneling conductance dl/dV at posi- 
tive and negative bias, the corresponding Bogoliubov an- 
gle <d{x,y), and the logarithm of the absolute value of 
its Fourier transform. We consider the following values 
for the bias: V = ±0At,±0.8t,±1.2t. Note, that at bias 
V = ±0.4i, that is under the gap value A « 0.8t, the 
pattern of the local Bogoliubov angle (see Fig. [6h is ro- 
tated 45 degrees with respect to the patterns calculated 
at higher biases V = ±0.8t (near the gap, see Fig. [T]:, 
and V = ±1.2t (above the gap, see Fig. [SJ:). The sites 
on the lattice where there is a large particle-like com- 
ponent of the Bogoliubov excitation, hole component is 
small. Complementary pattern is observed on opposite 
bias. This "rotation" is commonly present in the whole 
field of view. 

We also present BA along the diagonal line cut for 
our numerical calculation to compare with experimental 
results. We observe an out-of-phase angle change for low 
energy vs high energy BA, Fig. [HI This out of phase 
behavior is consistent with the behavior seen in Figd) 



VI. CONCLUSION 

In conclusion, we have introduced a new spectroscopic 
measure, Bogoliubov angle <d{ri,E). This measure al- 
lows one to image local particle-hole admixture in the 
superconducting state and in the normal state with su- 
perconducting correlations. 

Bogoliubov angle can be studied as a function of po- 
sition. It can also contain nontrivial Fourier compo- 
nents. This could allow us to make connection with the 
spatial interference of quasiparticles in superconducting 
statei^i^'^i^ Complementary to the momentum space 
information one can look at the energy dependence of 
BA. Energy dependence observed experimentally clearly 
indicates that there is a change in behavior in Q{r,E) 
at ~ 20 — 25mV, Figs l4l5l This energy range clearly 
correlates with the changes in the interference patterns. 
We interpret these changes as an evidence for a change 
in the superconducting coherence that is weakened at 
higher energies. 

As a future application we propose that BA be studied 
as a function of doping and temperature. One can use to 
investigate BA and particle-hole at temperatures above 
Tc for studies of the nature of the pseudogap phase. Us- 
ing Bogoliubov angle one could be able to identify how 
robust the particle-hole mixture is in the normal state 
and therefore be able to differentiate between different 
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scenarios of PG state. 
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APPENDIX A: Anderson Mapping 

Here we recall the Anderson— mapping of reduced 
BCS model on the effective spin model. The reduced 
BCS Hamiltonian is taken to be 



Complete identification with Bogoliubov quasiparticles is 
clear if one identifies 



sin 9k ~ 2wkWk, cos 6k = Wk ~ '^11 (^"8) 
To make a contact with BA we notice that in case of 
broken translational symmetry we can work out exactly 
the same representation based on eigenfunctions in real 
space u„(ri), z;„(ri). Then the mapping on the spin prob- 
lem will be done in real space, angle Ok is proportional 
to the BA defined in the Introduction, and we will have 
0£;(ri) as defined in Eq.([21). One can immediately see 
the direct connection with the Anderson angle used in 
this effective spin model. 



Hred = -X!^^'' ~ ^^^^ n_k) - 

k 

^k.k'C^clkC-k'Ck' = (A-1) 

kT^k' 

-2 ^(ek - ^i)s.M - 1/2 5I(4*k' + (A-2) 

k k,k' 

where we assumed translational invariance for simplicity 
and omit spin indexes. Spin operators are defined as: 

Sz,k = 1 - "k - "ks (A-3) 
4 = K=^UU, (A-4) 
Sk = ^k = CkC_k (A-5) 

and they represent a complete spin algebra over space 
rik — k = 0, the so called hard core boson constraint, 
z component of the spin corresponds to state with well 
defined particle number and corresponds to pairing 
correlations. Anderson showed that this reduced Hamil- 
tonian describes the spin Sk in an "external" field point- 
ing at an angle Ok 

e, = 1/2^-'^--'""^- (A-6) 
Ck - M 

One immediately recognizes this a self-consistency equa- 
tion for BCS solution, once we assume Vkk' to be constant 
in a range near Fermi surface. Excitation spectrum for 
the effective spin model is: 

Ei, - ((ek - m)' + 1/4(5] lik' sin Ok')')'/' (A-7) 

k' 



Appendix B: Numerical details 

The numerical solution of Eq. ipO)) together with the 
self-consistency condition Eq. ([2T|) requires iterative solu- 
tion and it is organized as follows: 

(1) For a reasonable initial value of the order parameter 
As{ri) we solve the eigenproblem Eg. ipn)) to obtain the 
quasiparticle amplitudes (M„(ri), u„(ri)) and the quasi- 
particle spectrum En- 

(2) Then, substituting (u„(ri), w„(ri)) and En into 
Eq. (|2ip we compute a new approximation of the order 
parameter A^"^^'^^ (r^). 

(3) In order to avoid numerical instabilities during 
iterations, we use a mixing scheme A^""*"^' (r^) = 

aA^°^^''^(r,) + (1 - a)A^"'(r,), where A^"^(r,) is the or- 
der parameter at the previous iteration step. Adjustable 
parameter a is a number between and 1. To insure 
convergence, we increase the current value of a by 5% 
if the relative deviation between two consequent steps 
5" = max,,5 |A^(r,)- A^-i(r,)|/max,^5 |A^(rO| has de- 
creased, < S*"^^. And we decrease a by 20% in the 
opposite case, S"' > S^^^. 

(4) The computed A^"^^-'(ri) is used for the next itera- 
tion step. 

We repeat iterations until we achieve the acceptable 
level of accuracy (e = lO^''). After the end of the pro- 
cedure, we perform an additional step with a = 1 to 
ensure convergence of the obtained solution. It usually 
takes 20 — 40 iterations to converge. 
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FIG. 2: (Color online) BCS coherence factors u'^(k) (blue), 
w^(k) (red) are shown as functions of energy. The function 
C(k) = |u^(k) - i;^(k)| = | cos2e(k)| (black), shows substan- 
tial departures from unity only in the energy range on the 
scale of the gap A near the Fermi energy, where there are 
substantial pairing correlations. 
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FIG. 3: (Color online) g{'r,V) with 54nm FOV at V =-18mV (a) and V =+18mV (b) with their fourier transforms in the 
insets. The modulations are visible and consist of several wave vectors, c, Typical averaged spectra taken at different area, d, 
same spectra as c, but zoomed at the low energy feature below the maximum gap. Systematic deviation in spectra between 
the negative and the positive sample bias is seen, as indicated by arrows. 
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FIG. 4: (Color online) a, 6(ri, V) with 54nm FOV at V =18mV. b, fourier transform of 6(ri, V) in a. c, Distribution of 9(ri, V) 
at V =18mV. d, Spatial evolution of the fourier filtered dl/dV at 18mV and -ISmV (black, darker) and 0(ri, V = 18mV) (red, 
lighter) with 2iT/qr modulation along the red (light solid) line starting from the solid (red) circle in a. 
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FIG. 5: (Color online) a, b, c, d, e, images of the BA at each bias voltages (V^ =10, 18, 26, 34, and 42mV) and their fourier 
transforms, f, distributions of the BA at each bias voltages from to 90mV. Peak positions of the histogram are traced by 
black line. 
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FIG. 6; (Color online) a,b,c,d, calculated LDOS on a square 
32 X 32 lattice at T = 0. We assume 40 randomly placed 
impurities with individual impurity strength V""p = it. The 
pairing strength is set Vint = —2t and chemical potential 
fj, = (see text), a, calculated local ^ tunneling conductance 
at positive bias V = OAt. b, calculated local ^ tunneling 
conductance at negative bias V = —OAt. c, corresponding 
Bogoliubov angle Q{x,y). d, logarithm of the absolute value 
of Fourier transform of BA logj^o(|0(fca;, ky)\) with subtracted 
average value < Q{x,y) >— 45°. 
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FIG. 7: (Color online) a,b,c,d, calculated LDOS on a square 
32 X 32 lattice at T = 0. We assume 40 randomly placed 
impurities with individual impurity strength V""p = it. The 
pairing strength is set Vint = —2t and chemical potential 
fj, = (see text), a, calculated local ^ tunneling conductance 
at positive bias V = 0.8t. b, calculated local ^ tunneling 
conductance at negative bias V = —0.8t. c, corresponding 
Bogoliubov angle Q{x,y). d, logarithm of the absolute value 
of Fourier transform of BA logj^o(|0(fca;, ky)\) with subtracted 
average value < Q{x,y) >— 45°. 
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FIG. 8: (Color online) a,b,c,d, calculated LDOS on a square 
32 X 32 lattice at T = 0. We assume 40 randomly placed im- 
purities with individual impurity strength = it. The 
pairing strength is set Vint = —2t and chemical potential 
fj, = (see text), a, calculated local ^ tunneling conductance 
at positive bias V = 1.2t. b, calculated local ^ tunneling 
conductance at negative bias V = — 1.2t. c, corresponding 
Bogoliubov angle Q{x,y). d, logarithm of the absolute value 
of Fourier transform of BA logj^Q{\Q{kx, ky)\) with subtracted 
average value < Q{x,y) >= 45°. The Fourier transform we 
obtain is consistent with FT intensity us seen in the exper- 
iment, see inset in Fig. [S]e. Please note rotation of {qx,qy) 
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FIG. 9: (Color online) Profile of the Bogoliubov angle O 
along the line cut for bias values 0.4t, red (dashed) line; and 
1.2t, black (solid) line. The line cut is taken along the di- 
rection [1,1]. Note the angle inversion effect with respect to 
the optimal mixing angle value of 45° for low energy and high 
energy BA. 



